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ABSTRACT 

We generalize Doroshkevich's celebrated formulae for the eigenvalues of the initial 
shear field associated with Gaussian statistics to the local non-Gaussian f n i model. 
This is possible because, to at least second order in distributions at fixed overden- 
sity are unchanged from the case f n i = 0. We use this generalization to estimate the 
effect of f n i ^ on the abundance of virializcd halos. Halo abundances are expected to 
be related to the probability that a certain quantity in the initial fluctuation field ex- 
ceeds a threshold value, and we study two choices for this variable: it can either be the 
sum of the eigenvalues of the initial deformation tensor (the initial overdensity) , or its 
smallest eigenvalue. The approach based on a critical overdensity yields results which 
are in excellent agreement with numerical measurements. We then use these same 
methods to develop approximations describing the sensitivity of void abundances on 
fni. While a positive f n i produces more extremely massive halos, it makes fewer ex- 
tremely large voids. Its effect thus is qualitatively different from a simple rescaling 
of the normalisation of the density fluctuation field as- Therefore, void abundances 
furnish complementary information to cluster abundances, and a joint comparison of 
both might provide interesting constraints on primordial non-Gaussianity. 

Key words: methods: analytical - large scale structure of the universe 



1 INTRODUCTION 

Detections of non-gaussianity can discriminate between different inflation models (e.g. Maldacena 20031 ). The local /„j model, 
where the primordial perturbation potential is 



* = 4>+/n^ 2 -^ 2 }), (1) 

where </> is a Gaussian poten t ial field and f„i is a scala r, has been the subject of much recent study (e.g., Buchbinder et ail 
2008; Khourv fc Piazzall2008l: ISilvestri fc Trodden 20081. and refer e nces therein). Constraints on this model tend to be of tw o 



types - from the CMB (jrComatsu et al.l 1 20091 : Ifiikaee et aljEooi: lYadav fc Wandeltlboosl: 



200 



^) and from large scale structures in the Universe ( Kovama*etldTfl999f: 



Matarrese et al 



McEwen et al 

IT 



2000n Scoccimarro et al 



20081: 



Rossi et al 



200. 



4; 



Sefusatti fc Komatsul2007l;llzumi fc SodalpOollLo Verde et al.ll2008l:lDalal et al.ll20q8l;lMatarrese fc Verdel2008l;ICarbone et all 

2008 : Afshordi fc Tollevll2008l;lsiosar et alj2008l:lMcDonaldll2008l: Taruva et al.ll2008l;lsiosar 2009l : lGrossi et all2008l ;l Kamionkowski et al 



20091 : IPillepich et alJliool IPesiacques et alJIiooi : lLam fc Shethll2009bl : iGrossi et al.ll2009l ). 

The initial shea r field is expected to play an important ro l e in the formation of la r ge scale structure s ( Zel'Dovichl 197o1 : 
Bond fc Mverslll996l : Eee fc Shandarinlll998l : Isheth et al-lEooH : iBernardeau et al.ll2002llDesiacauesll2008l ). The main goal of 
the present work is to show that much of the machinery developed for the study of structure formation from Gaussian 
initial conditions can be carried over, with minor modifications, to the study of /„; models. We do so by showing how to 
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generalize Doroshke vich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics 
l Doroshkevicblll970h to the local non-Gaussian /„j model. 

We then study how the abundance of virialized dark matter halos depends on / n j . This study focuses on two problems: one 
is the physical model for halo formation, and the other is the statistical pr oblem of how this collaps e model is used to estimate 



the ab undance of such collapsed objects. We use the approach pioneered bvlPress fc Schechterl |l97J) and refined bv lBond et al 
1 1991) to address the statistical problem (see Maggiore fc Riottoll2009a ; Lam fe Shethll2009a . for recent extentions which treat 
the f n i 7^ case). We study two different mode ls for the physics of halo formation: one in which halos form from sufficiently 
overdense regions in the initial fluctuation field ( Press fc Schechterlll974 ; Sheth et al. 2001 ) , and another in which the criterion 
for h alo formation is that all three eigenvalues of the ini tial deformation tenso r exceed a certain value (e.g. iLee fc Shandarin 
1998) . In the form er, we use the triaxial collapse model of Bond fc Mversl 1 19961 ) to estimate this critical overdensity (following 



Sheth et al.ll200ll ). If our way of estimating halo abundances are reliable, then comparison with simulations run for a range 



of fni values provides a novel way to study the physics of gravitational clustering. 

Sect ion [2] provides explicit expressions for the initial eigenvalue distribution, and for the initial distribution of the variables 
which arise naturally in the context of triaxial collapse models. These are used, in Section [3] to estimate how the mass 
function of virialized objects is modified when /„; 7^ 0. This Section also shows the re sult of comparing these estimates with 
measurements in simulations. With some care (e.g. ISheth fc van de Weyga ert 2004), much of the analysis can be carried 
over straightforwardly to study void abundances; this is the subject of Section 0] A final section summarizes our results. An 
Appendix describes an alternative estimate of halo abundances which is logically consistent with previous work, but which 
does not reproduce the /„; dependence seen in simulations. 



2 THE LOCAL NON-GAUSSIAN MODEL 

We are interested in models where the primordial perturbation potential is given by equation {l}. We will use P<t,(k) to 
represent the power spectrum of fa in what follows we will set P<t>(k) = Ak ns ~ 4 , where n s ~ 1, and A is a normalization 
constant that is fixed by requiring that the rms fluctuation in the associated non-Gaussian initial density field (which we will 
define shortly) has value eg. The power spectrum and bispectrum of the $ field are 

P 4> (fc) = P,(fc) + Jdq [P*(q)P+(\k -q\)- P*(k)P+(q) - P^tyP+Qk - q\)} , (2) 

B#(fci,&2,M = 2/ni [P^(fci)P«(fca) + cyclic] + 0(f* l ) (3) 



( Scoccimarro et al.l 2004) 



2.1 The shear or deformation tensor 

Define D as the real, symmetric 3x3 tensor whose components are proportional to the second order derivatives of the 
potential $: 

$ij = fa + 2f nl (fafa. + 4><t>ij), (4) 

where fa — dicj) and faj = didjfa We will sometimes refer to D as the shear field associated with the potential Correlations 
between the <I>.y will be very useful in what follows. These depend on the correlations between <f> and its derivatives but, 
because <f> is Gaussian, they can be computed easily. Doing so shows that the six components of D are not independent: 
although the three off-diagonal terms are not correlated with the others, the three diagonal terms are. However, if we set 

x = ^2$ii, y= i ($ lx -$22), z = £ ($11 +$22 - 2*33), (5) 



then these three parameters, combined with ($12, ( I > 23, $31), form a new set of six independent components (|Bardeen et al 
19861 ). When / n ; = 0, then each of these is an independent Gaussian random field. 



Most of the complication in /„; models arises from the fact that we are almost always interested in spatially smoothed 
quantities. Fortunately, smoothing is a linear operation, and the new variables x,y,z axe just linear combinations of the 
elements of D. Hence, if W(kR) denotes the Fourier transform of the smoothing window of scale R, to second order in /„;, 

(x 2 ) = ° 2 , {y2) = i> {z2) = i' <*«f>*j = ^ (6) 

(x 3 ) = 2/ n!7 3 <y 3 > = </>=0 <$ 3 J ) l#J = 0, (7) 
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where 



(2tt) 3 
2 



dk 
~~k 
dfei 
"fa" 



4tt k 7 M 2 (k) P<s,{k) W 2 (kR), 
dk 2l 



klM(ki)W(kiR) 



-/ 



k a 2 M{k 2 )W{k 2 R) / d/xia fel 2 M(&i 2 )W(fci 2j R) 



B i ,(k 1 ,k 2 , k 12 ) 

2fnl 



(8) 
(9) 



and M(k) = (3D(,z)c 2 )/(5f2 m -ffo ) T(k), where is the CDM transfer function and -D(z) is the linear growth function. In 

what follows, the quantity aS'3 = a (x 3 ) / (x 2 ) 2 , will play an important role, bec ause it represe n ts the l eading order contribution 
to the non-Gaussianity (note that it is proportional to /«;). Appendix A in lLam fe Sheth (2009a) provides a useful fitting 
formula for this quantity. 



2.2 Joint distribution of eigenvalues 

Equation shows that, to first order in /„;, five of the six parameters have zero skewness, so, to first order in /„;, all but x 
are drawn from Gaussian distributions. However, x is the trace of D, so it is the linear theory overdensity S. If p(S\R) denotes 
the distribution of 8 when smoothed on scale R, then the fact that the other variables have the same distribution as in the 
case fni = allows one to provide an excellent analytic approximation to the joint distribution of the eigenvalues. Namely, 



p(X\R) = p(5\R) j exp + — j (Ai - A 2 )(A 2 - A 3 )(A 1 - As), (10) 

where 5 = Ai + A2 + A3, I = A1A2 + A1A3 + A2A3, and recall that a is a function of R. Our convention is Ai > A2 > A3. This 
has the same form as Doroshkevich's (1970) formula for Gaussian fields; the only difference is that here p(S\R) is not Gaussian 
(we provide an expression for it in equation 1131 below) . 

The fundamental reason why this works is that Doroshkevich's formula is actually the product of two independent 
distributions, one of 5, and the other of a qua ntity which is a combination of the five other independent elements of the 



deformation tensor (e.g. ISheth fc Torm cn 2002). Since the distribution of each of these other elements is unchanged from 
the Gaussian case (we just showed that they all have zero skewness) , this second distribution is unchanged from that of the 
Gaussian case - the only change is the dist ribution of 5. In fact , this h olds for any local mapping $ = f((j>) (and is also true 



for the alignment of the principal axes, see Desiacques fc Smith ( 20081 )) 



2.3 Distributions at fixed S 

One consequence of this is that distributions at fixed 5 are the same as in the Gaussian case. For example, 

P ^ X ^ = rl^y (2^) exp ("S + 1^) (Al " Aj)(Ai + 2X3 - S){2X * + Aj ~ S) = P°( Xi ' x ^' ( n ) 

where i 7^ j can take values from 1 to 3, Uj = A^Aj + (Ai + Aj)(<5 — Ai — Aj) = S (Ai + Xj) — (\ 2 + AiAj + X 2 ), and the 
subscript indicates the distribution associated with Gaussian initial conditions, for which /„j = 0. Integrating over one of 
the eigenvalues in the expression above, e.g., Ai, yields an expression for the distribution of the other at fixed 5. Clearly, such 
expressions will also be the same as in the Gaussian case. Hence, 

p(A#)= P0 (A#) (12) 

far j = 1,2, 3. 



2.4 Edgeworth approximation for p(S\R) 



Because we are interested in small departures from Gaussianity, the Edgeworth expansion provides a convenient form for the 
distribution of 8: 



p(S\R) d6; 



1 ( a(R)S 3 (R) H ^ 



6 



where a(R) is given by eq uation (El), a S 3 



a(R) 

{x^/ix 2 }^ 2 



■ d5 = 



2na(R) 

2/ ni7 3 /^ 3 , and H a (v) 



1 , TT I 

1 + —H z (v 



Po(S\R) dS, 
v(y 2 — 3) with v 



(13) 



S/a(R). (See 

Lo Verde et al.l l|2008l ) and lLam fc Shethl l|2009bl ) for previous work with the Edgeworth expansion in the context of / n ; 



models.) The final equality writes p as a correction factor times the Gaussian distribution po to highlight the fact that 



p(X\R) = 



1 + (S/a) 



P0 (X\R) 



(14) 



where po(X\R) is Doroshkevich's formula. 
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2.5 Distribution of eigenvalues 



Replacing p(5\R) in eg . (110(1 by its Edgeworth expansion and integrating over two of the three eigenvalues, we can write 



p(Ai) 



= Po(Ai) + ^Ap(Ai) 



where 
Po(Ai) 



Po(A2) 
Po(A3) 



and 
Ap(Ai 



Ap(A 2 ) 
Ap(A 3 ) 



Yl-KG 



Ai / 9A?\ 



20 — exp 
a 



2tt exp 



+3v37rexp 



2 v / ¥ct 
V5 



exp 



4(T 2 



127TCT 



20— exp 



15Af 

4cr 2 



9A§\ 
2a 2 I 



1 +erf 



5Af \ 
2a 2 J 

v^Ai 
2a 



X 2 

1-20^| 
a* 



3v37rexp 



4a 2 



27r exp 
v^Aa 



|*|WcfV2^H 1-20^| 



2a 2 



erfc 



2a 



VE 

YI-ko 



+3V3tt 



435Aj ! 
4a 2 

25Ai 
4a 



9A? 
2a 



25Ai 
21a 



\2 i4 

51 - 185^r + 100-^ 1 exp 



5A? 

2(7 2 



— - — 1 exp 



5X3/2 



15A; 



(I 



V2a V 2 ^ 2 



3 1 po(A2) = 



15A 2 
4a 2 

3/2 



#3 



erf 
15A 2 



\/2a 



( V3Xi 
V 2a 

Po(A2 



127ra 



25 
27 



+3V37T 



435 A| Xi 
8-- 1 + 100^ 

4 cr 2 cr 4 

25A 3 /5A| 



exp 



4cr 



2a 2 



exp 



2a 2 
15A| 

4(T 2 



25A 3 



erfc 



27a 

/ \/3A 3 
V 2a 



^3 , „„ A| 



51 - 185-^ + 100 



5Af 
2a 2 



exp 



^3 

2a 2 



(15) 



(16) 



(17) 

1+crf f >/2 — 
a 

(18) 
(19) 

erfc f a/2- 



(20) 



These results should be easily extended to the distribution of shear eigenvalues at multiple points (e.g., Desiacqucs & Smith 
I III)). 



2.6 Distribution of 5, e and p 



In addition to the individual probability distriburions of the three eigenvalues, we can also derive expressions for the quantities 

of most interest in the ellipsoidal collapse model. These are the ellipticity e and prolateness p, where 

Ai — A3 . Ai + A3 — 2A2 1 3A2 A2 — A 3 

_ =e 



e = 



and 



25 ' ~™ P 25 2 25 5 

The discussion above means that the distribution of e and p at fixed 5 are unchanged from the Gaussian case 

5 



1 1125 , 2 

g(e,p\5) = 7== e(e 



' 10tt 



P 



However, the distribution of ellipticity is changed: 



9(e) 



dp / d<5c?(e,p|o)p(o) 



-(5/2)(J/<r)- , C3e- 5 +P'') 



,70(6) + ^A 5 ( e ), 
6 



where 



and 



9o(e) = 



45e 



(l + 20e 2 )(l + 15e 2 ) 5 / 2 



VEe(l + 30e 2 )\/l + 15e 2 - (1 + 20e 2 ) arctan 



VEe 



VI + 15e : 



-45000 e 4 

, |-r e 

vTutt e 



4725e b + 90e 4 - 26e 2 - 1 



(21) 



(22) 



(23) 



(24) 



(25) 



(l + 15e 2 ) 4 (l + 20e 2 ) 5 /2 

To check these expressions, we have used a Monte Carlo method to generate (x, y, z, $12, $23, $31) (all but x are drawn 
from Gaussian distributions). We then solve the eigenvalue problem (by solving the associated cubic equation) to obtain 
(Ai, A2, A3) and hence (5i, e,p). We then compute the distributions of the eigenvalues, and of 5, e and p, and compare them 
with the associated quantities when /„; = 0. The symbols in Figure [T] show our Monte Carlo results when the smoothing scale 
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Figure 1. Difference between initial distributions in the /„; and Gaussian models for f n i = 100 (lower panels) and /„; = —100 (upper 
panels) when the smoothing scale is lfc -1 Mpc. In both panels, the right and left plots show p(Ai/cr) (orange curve, triangle symbols) 
and p(\^/cr) (green curve, hexagonal symbols), respectively. For clarity, these curves have been shifted to the right and to the left by 3. 
The middle plot shows P(X'2/t) (blue curve, square symbols). 



is l/i _1 Mpc, and the smooth curves show the analytic formulae derived above. Notice the reflection symmetries in Figures [T] 
and [5] of the opposite sign of /„;: it is due to the fact switching the sign of /„; only changes the sign of aSz/Q without 
modifying other terms in Ap(Xi) and Ag(e). 



3 HALO ABUNDANCES 



Press &i Schechterl (|l974f ) (hereafter PS) argued that the abundance of collapsed virialized halos may be estimated from the 



statistics of the initial fluctuation field. They used the assumption that halos form from a spherical collapse to argue that 
such objects started out as sufficiently overdense regions in the initial fluctuation field. They then used Gaussian statistics to 
estimate collapsed halo abundances. Although the way in which they used Gaussian statistics to make the estimate is flawed, 
Lee fc Shandarin (1998) (hereafter LS) suggested that one might be able to provide a better estimate of the abundance of 



collapsed halos by repeating the PS argument, but changing the collapse model to allow halo formation to be nonspherical. 
In particular, they suggested that one should identify halos with regions in the initial field where all three eigenvalues were 
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Figure 2. Difference between initial distributions of ellipticity e in the /„; and Gaussian models for /„; = 100 (lower panels) and 
f n l = —100 (upper panels) when the smoothing scale is l/i -1 Mpc. 

greater than some critical value, A c . We will use the analysis above to show how this estimate of the halo mass function 
depends on /„;. 



3.1 If halo formation depends on the initial overdensity exceeding a critical value 

The PS-like estimate of the mass fraction in halos above mass M is 

roc roo V ci 

F(> M) = F{< a{R)) = J dSp(S\R) = J dS 1 + -jptf 3 Po(5\R), (26) 

where R = (3M/47rp) 1 ^ 3 , so a(R) is actually a function of M, and S c is the critical density required for collapse in the spherical 
model. So, 



and hence 



dF 8F d(aS 3 /G) , _i» , s aS A dhiv c 

■5- = -5 r- 5 (y c - v c )v c p {v c ) v c p (v c )H 3 (v c ), 

da da da 6 da 

OF dF ( aSz „ , , , -iv d(aS s /6) \ dF ( ,crS 3 . 
d In v c a vsxVc \ 6 omVc / dmu c \ 6 



(27) 
(28) 
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Lo Verde et al 



2008); the term in brackets is the ratio of the halo mass function when / n ; 7^ to that when / n ; 



Lam fc Shethl l2009al . argue that, formally, the term in brackets is not the full story, but that it is, nevertheless, a good 



approximation.) The ratio (OF /din i/ c ) /(dFp/d ln u c ) provides a good description of the fractional change in the halo mass 
function induced by the coupling parameter /„; ( Desiacaues et al. 20081 ) . even though the functional form of dFo/dln v c does 
not provide a good description of halo abundances in Gaussian cosmologies. 

Halo abundances in /„; = simulations are usually well-approximated by the functional form of Sheth fc Tormen ( 19991 ): 

apST - aI / 2 /2 
v f(y)=^^=2A^v [l+(^z/)- 2p ]^^, (29) 



d In v ' L v " V27T 

where p = 0.3, a — 0.7 and A = 0.322 comes from requiring that the integral over all v equal un i ty. Eq uation (|26|l with /„; = 
would yield a = 1 and p — 0, so A would be modified appropriately. Recently, iGrossi et al.l (J2009|) have studied the effect 
on the f n i model of s imply setting 5 C — > \/a8 c in equation (|26p . In effect, they ignore the consequences of p 7^ 0. However, 



Lam fc Shethl (|2009al ) have shown that when the non-Gaussianity is weak, then the correction factor for halo abundances is 
well approximated by 



dF/d\nu c 
dFo jd In v c 



6 



B(a) 



where 



(30) 



(31) 



with a — 0.7, ft = 0.485, and a = 0.615, is motivated by models of triaxial collapse (|Sheth et al.ll200ll ; ISheth fc Tormenll2002l ) . 



3.2 If halo formation depends on all three initial eigenvalues exceeding a critical value 



Suppose that the criterion for halo formation is not that S, the sum of the eigenvalues, exceeds 5 C , but that the smallest 
eigenvalue A3 exceeds A c . Then, the analogous argument yields 



5/3 

F(> M) = F(< u(R)) — I dSp(S\R) I d\ 3 p(\ 3 \5)= I dS 



roc t* 

/ d5p{5\R) / 



1 + 

6 



where 
Po(A 3 > \c\5) 



3V10 (S - 3A C ) 



exp 



5(5 -3A C ) 2 

8(7 2 



4 v / 7r cr 
If we define 4 = \ c /a(R), then 

F(>4) = *o(> 4) + AF(> 4), where 

6 



+ - 2 1 erf 



(5-3A c )V10 



4cj 



+ erf 



p (S\R)P (X 3 > \ C \S), (32) 



2(7 



(34) 



AF(> 



d^ exp( -Z 2/2) ^ 2 -3)Po(.-34) 



3la 

5 3/2 



2tt 



162V27T 



(100£e - 105^ 2 + 9) exp(-5^ 2 /2) erfc(V24) 



-3/2 



48^ 



(2 - 15£ 2 C ) exp(-15^/4) erfc( v / 34/2) 



+ig^4 exp(-9£ 2 /2) (5-8£ 2 ). 
6487T 



(35) 



The halo mass function is 
OF 



8Fo d(aS 3 /6) f) ,*S 3 d(AF) 

54 + d£ c lc) + ~e dlV 



dF 

d£ c 

d(AF) 



no (u\ 1 

25 VE 



2 2 4 3 4 tt 



3 12 

exp 



exp 
94 



5% 
2 



erfc(v / 24 



4^ 



exp 



154 



erfc 



^/3£ c \ ->\ r < 



+ 



(64 - 8704 + 8004) - V27r84(51 - 1854 + 1004) exp 



3tt 

54 



£ c exp I — 



94 



erfc(v / 24) + 



/ 3^3 4 4(2-54) exp 



erfc 



^ n/34 



Notice that 



d£ dFo 

d£ c * 94 



( 1 + ^^4)-(4-4-)») 



(36) 



(37) 



Thus, in this case, the /„; modification to the halo mass function is qualitatively different from that associated with the 
spherical evolution model. 

We mentioned above that the ratio (dF/dhi v c )/(dFa/d\n u c ) is well described by the term in brackets in equation (|28[) . 
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Figure 3. Halo multiplicity function dF/d\nv for the Gaussian simulations at z = 0.3. Black symbols show the measurements from 
simulations; short-long-dashed, dotted, dot-dashed and solid curves show equations {26}, 1136 ft . 1291 . and 129^ with a new normalization 
(A = 0.26) respectively. 



Since v c and la are linearly proportional to one-another, it is interesting to ask if (dF/dln£ c )/(dFo/d\n£ c ) is also well 
described by the term in brackets in equation (|28p . (We have already shown that this ratio is not described by simply 
replacing v c — > 3£ c -) Figure [4] shows that it is not. Thus, if we were certain that the logic which leads to this estimate of the 
mass function were reliable, then we would conclude that, by studying how halo abundances depend on /„;, we may have 
learnt something important about the physics of halo formation: the initial overdensity matters more than the value of the 
smallest eigenvalue. 



3.3 Comparison with measurements from numerical simulations 

Figure |3] shows the multiplicity function measured in the f„i = simulations of Desiacaues et all (20081 ). where a detailed 
description of the runs can be found. (Our analysis is complementary to that of ( Pillepich et alj 20081 ). who have recently run 
a large set of simulations of the f n i model; they studied halo abundances and clustering in their simulations.) Curves show 
equation (|26l) (magenta, short-long-dashed), equation ()36[) (black, dotted), and equation (|29|) (cyan, dashed) with S c = 1.66 
and Ac = 0.41. In contrast to most previous work, equation (13611 appears to g ive a better fit than equation (I29[) : this may 



be due to the fact t hat the halo finder (AHF, see iKnollmann fc Knebe (2009)) is not standard. (The halo-finder used by 
1 200ct ) is more standard, and they indeed find that equation 1291 with standard choices for its free parameters, 
works well.) The red solid curve shows equation ([29} with a new normalization (A = 0.26) and the agreement with the 
numerical measurements is much better. 
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0.6 0.8 1 2 4 6 0.6 0.8 1 2 4 6 

v = 5 / a v = 6 / a 



Figure 4. The ratio of halo mass function of f n i = 100 models to the corresponding Gaussian models using equations ((28} , (30jl, iTAl")) 
and (1361 . (magenta short-long-dashed, solid red, cyan dot-dashed, and black dotted). Panel on left shows /„; = 100 and panel on right 
shows f n i = —100. 



It is conventional to show the effects of f n i 7^ on the mass function by ratioing with respect to the /„; = case. The 
symbols in Figures [4] show this ratio for /„; = 100 (left) and —100 (right) respectively. The short-long-dashed, dotted and 
dot-dashed curves show these ratios for the same models as in the previous figure (equations 1281 1361 and IA1|I . The solid (red) 
curve shows equation (|30[) . The est imate motivated by the spherical collapse model (equations | 28l ) describes the measured 
ratio very well (in agreement with Lo Verde et al. 20081 ; Desiacques et al. 20081 ; Grossi et al. 20091 ). even though the mass 
function on which it is based is a bad fit to the /„; = data. On the other hand, the same logic applied to a prescription 
based on the smallest eigenvalue (equation I36[) fits the Gaussian (i.e. f„i = 0) mass function reasonably well, but do es not 
describe deviations from non-Gaussianity very well! In addition, the same logic applied to the Sheth fc Tormen ( 19991 ) mass 
function (equation IA1|) also does not fit the ratio very well. However, our excursion set based approach (equation I30|l seems 
to match the measurement as well as, if not better than, any of the other methods. Note that it is the only model which 
mat ches both the f n i = mass function, and the /„; 7^ ratio. We also checked that our excursion set based approach agrees 
with Pillepich et al.l 1 20081 ) fitting formula's prediction (less than 6% difference in the range of validity of the fitting formula, 
which is 1.4 < v < 5). 



4 VOID ABUNDANCES 

Underdense regions are also a good probe of / n ; ( Lam fc Sheth 2009b ). Kamionkowski et al. ( 20091 ) have applied the analog 



of equation (|26[) to study void abundances when /„( 7^ 0. In what follows, we estimate void abundances associated with the 
analog of Section 13.21 b y assuming they are simply the opposite of halos. However, because of the 'void-in-cloud' problem 
identified bv lSheth fc van de Weygaertl (J200J), even when /„; = 0, these analyses are, at best, appropriate only for the largest 
voids. Hence, we also discuss the effect of including the correction for the 'void-in-cloud' problem. 



4.1 If voids formation depends on all three eigenvalues lying below a critical value 



Given a critical value A„, the corresponding mass function of voids is 

dSp(8\R) / dA 1 p(Ai|<5) = / ,M 

-00 J5/3 

where 

5(3A„ -<5) 2 ~ 



1+^3 
D 



, 3V10 (3A„ - S) 

Po(Xi < X V \S) = <^ - exp 

1 4y7T a 



8a 2 



+ 5 I-' 



(3A„ - S)VW 



4(T 



Po(5\R)P (Xi < X v \5), 
(3X V -5)VlO 



erf 



2(7 



(38) 



e(3A u -5). (39) 



10 T. Y. Lam, R. K. Sheth & V. Desjacques 



If we define £ v = X v /a, then the corresponding Fq(< £ v ) and AF(< £ v ) can be obtained by replacing erfc(s) by —[1 + erf(a;) 
in equations (|34fl and (|35[) . The void mass function can be obtained by making the same replacement, so it is 

aS 3 d(AF) 



OF _8F d(aS 3 /6) 
dl v dl v d£ v { - 



dF 



10 {51. 



djAF) _25 y/E 

d£ v 



2 2 4 3 4 tt 



1 

~ 12 
exp 



exp 

%£ v 
~~2 



») + ■ 

5£? 
2 



6 



») [l + erf(^)]+^g 
J (64 - 870£ 2 V + 800^) + v^tt84(51 - 185^ + 100^) exp 
1 + erf 



: exp 



15lg 
4 



1+erf 



37r3 4 ^(2-5^)exp 



( V3£y 

V 2 



, 5^5. 
+ -= — 4 exp 
Sn 



Q£ v 
~2 



5£ v 



l+erf(V24,)l 



15£g 
4 



^ 734 



(40) 



4.2 Void-in-cloud problem: excursion set approach 

An important aspect in the void abundance is the overcounting of the voids located inside collapsi ng regions. The formalism 
of counting voids as regions belo w some critical value (den o ted S v in iKamionkowski et al, (|2009h and A„ in the discussion 
above) does not account for this. ISheth fc van de Weveaertl (j2004h examined this problem using the excursion set approach 
by studying a two barriers problem: <5 C for haloes and 8 V for voids. 

We will now extend the calculation of the void-in-cloud effect to models where f„i 7^ 0. We use the constant barriers 
to demonstrate the method. Denote the two constant barriers correspond to the formation of halos and voids by 8 C and 8 V , 
F(s,S v ,5 c ) as the probability of a random walk crossing the barrier 5 V at scale s and it did not cross the other barrier S c . 
This probability is directly connected to the void abundances, including the void-in-cloud effect, as the random walk never 
crossed the halo formation barrier. It is related to the first crossing distribution f(s,S v ) by 

F(s,5 v ,S c ) =f(s,5 v )- [ S dSiF(Si,5 c ,5 v )f(s,5 v \Si,5 c ), (41) 
Jo 

where the second term on the right hand side substracts from the first crossing those walks that crossed 5 C at Si before 
crossing S v at s (but never had crossed S v before Si). Swaping S v and 5 C : 

F(Si,5 c ,5 v ) = f{Si,6 c )- dS2J r (S 2 ,S v ,S c )f(S 1 ,S c \S2,5 v ). (42) 
Jo 

Substituting equation (I42|) into equation (|4ip . we find 

F(a,6 v ,6 a ) = f(8,8 v )- f dSif(s,S v \Si,8 c )f{Si,5 )+ [ dSi dS 2 f(s, S v \Si, S c )f(Si, S C \S 2 , 8 V )F(S 2 , S v , S c ) (43) 



= f(.a,6 v ) + Y](-l) n / dSi... / dS n FT f(S m ,S m \S m +i,5 m+1 )f(S n ,5 n ), 

n=l Jo Jo m=0 

where the last expression is obtained after inserting equations (|41[1 and (|42[1 successively. Furthermore, So = s and 



(44) 



, 5 V if n is even 

° n = i x ■( ■ aa ( 45 ) 
o c it n is odd 



The nth order term in the summation of equation (|44|l denotes walks that have crossed the two barriers alternatively n times 
before crossing 8 V at s. Below we will work out the predictions of equation (|44[l for primordial Gaussian perturbations and 
for models with primordial non-Gaussianity of the f n i type. 



4-2.1 Gaussian initial conditions 



We would like to estimate f(s, 5 V \S, 8 C ), which is the first crossing probability of S v at scale s given it crossed the barrier S c at 
some scale S(< s). For Gaussian distributions with sharp-fc space filters, fo(s, 8 V \S, 8 C ) = fo(s — S,S V — 8 C ), and equation (|44l) 
reduces to 

~n-l 

Y\ fo(S m , Sm\S m + l,S m + l) 



Fo(s,8 v ,8 c ) = f (a,S v ) + J2(-l) n \ dSi 

n 

» fo(s, S v ) exp 



(•So 




/ dSi . . 


• / dS n 


Jo 


Jo 



fo(S n , S n ) 



5 C <iv 2 v 4 - 



(46) 



where the last expression is the approximation given by ISheth fc van de Weveaertl (|2004f ) with T> = \S V \/(8 C + \5 V \) and 
v = 5 v /s/s. 
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4-2.2 Local non- Gaussian f n i models 



The calculat ion of the condition al first crossing probability f(8 v ,s\8 c ,S) for the /„; model is analogous to that of halo 
abundances l|Lam fc She th 2009a). First we write down the probability p(s,8\S,8 v ) as 



provided that 8 < S v . So, 



P(s,S v \S, 5 C ) 



p{s,5\S,8 c ) = [ S dS' f{S',5 v \S,5 c ) P (s,5\S',5 v \S,S c ), 
J s 



dSp(s,S\S,8 c ) =/ dS' f(S',S v \S,8 c ) dSp{s,S\S',S v \S,8 c ) 



The derivative with respect to s is 

dP(s,S v \S, S c ) = f(s,6 v \S,5 c ) 
ds 2 



+ / dS' f(S\S v \S,8 c ) 



d /■»» 



ds 



&8p{s,8\S',8 v \S,5 c ). 



(47) 



(48) 



(49) 



The above equation is an integral equation for f(S v , s\S c , 5) and its zeroth-order solution is given by the left hand side of the 
equation (which can be evaluated using the bivariate Edgeworth expansion). We argue that, in analogy to the calculation of 
the halo abundance, the first-order solution is negligible compared to the zeroth-order. Therefore, we can make the following 
approximation 

,dP{s,S v \S,So) 



f(s,5 v \S,5 ) 



ds 



Note that for Gaussian distributions, 



dP (s,8 v \S,8 c ) d /-(*.-*e)A/5=S e"* 2 / 2 
— ' da; ■ 



ds 



ds 



2- 



-{8 V - 5 C ) exp[-(5„ - S c ) 2 /2(s - 5)] 



(50) 



(51) 



2(s-S) ^2tv(s - 5) 

which is the expected conditional distribution for Gaussian statistics. 

In Appendix[B] the right hand side of equation ((50} is evaluated using the Edgeworth expansion. Here, we will approximate 
the conditional first crossing probability by 

0-S3 



f(s,5 v \S,8 c )*s fo(s,6 v \S,6 c ) 



where 



(52) 



(53) 



ds \S V -8 C \ -^'"'v^ - 5 
Substituting equation (|52|l into equation (144 II and keeping only terms linear in (0-53/6), we can recast the integral equation 
into the form 

o53 



T(s,S v ,8 c ) = J r {s,5 v , 8c) + fo(s,8 v )H 3 — 



+ 



n=l J 



S„- 



dS n 



Y\ .fo(Sm,8m\S m +l,8m+l / 



fo(S„, 8 n ) 



y C(S m , 8 m , Sm+1, S m +l) + Hz 



'Sri 



(54) 



4.3 Comparison of models 

Figure [5] shows theoretical expectations of void abundances when /„; = 0: The solid (black) curve shows equation (|46ll 
with (8 V ,S C ) = (—2.81,1.66); The Long-dashed (cyan), dot-dashed (magenta), and short-long-dashed (green) curves show 
equation (|44fl . keeping terms up to n — 2,3,4 respectively. The zeroth order solution of equation (|44[) without the void-in- 
cloud effect is shown as the dotted (red) curve. It is the first crossing probability of a constant barrier and is the same as the 
PS prediction. Finally, the short-dashed (blue) shows the prediction of the LS formalism (equation I40p with A„ = —0.69. 

The predictions from the PS formalism and the LS formalism (both without the void-in-cloud effect) are different over 
a large range of void size, so comparisons with numerical measurements could distinguish which models describe best the 
formation of voids. Furthermore, while the effect of void-in-cloud is significant for small voids, it is negligible for big voids. 

The approximation equation (|46p provides a very good description of the solution to the integral equation (|44|l . even in 
the regime of very small voids (u 2 ~ 0.5) for which equation (|44[) requires the inclusion of high order terms. If one characterises 
the accuracy of the n-order term by the smallest v at which the inclusion of the next order term modifies the result by less 
than 1%, then the second and third orders are accurate for v > 0.9 and v > 0.6. 

Figure [6] compares the ratio of the void abundances for f n \ 7^ relative to the case f n \ = for the various analytic 
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Figure 5. Predicted void abundances for Gaussian initial conditions. Solid (black) curve shows equation 1 146 I I with (8 V , S c ) = (—2.81, 1.66). 
Long-dashed (cyan), dot-dashed (magenta), and short- long-dashed (green) curves are the solution to equation (144 It with n = 2,3,4 
respectively. Predictions without the void-in-cloud effect are also shown: dotted (red) is the zeroth order solution of equation 144L which 
is equivalent to the PS result; short-dashed (blue) shows the prediction from the LS formalism (equation l|40|l with X v = —0.69). 



approximations described above. The curve labels are identical to those in the previous figures, except for the solid (black) 
curve showing equation Q280 an d an addition dot-long-dashed (orange) curve showing the same equation upon neglecting the 
d(aSg)/ds term. These two cur ves, as well as the sho rt-dashed (blue) curve, include the scale dependence of it S3 using the 
approximation formula given in lLam fe Shethl (12009ar i. The other curves (dotted (red), long-dashed (cyan), dot-short-dashed 
(magenta), and short-long-dashed (green) for n = 0,2,3,4 in equation I44f) assume a constant a S3 = 0.022 As we can see, 
unlike halo abundances, a positive value of f„i increases the relative number of big voids, whereas a negative /„; decreases it. 

The overlapping of the solid (black), dot-long-dashed (orange), and dotted (red) curves justifies a posteriori our assumption 
of constant CTS3. For big voids, the dotted curve deviates only slightly from the other two where it agrees reasonably well with 
higher order solulions to the integral equation (|44[) . Hence, the difference between the various curves (apart from the blue 
short-dashed) are mostly governed by aS3 rather than the different functional form of the models. Note, however, that the 
prediction of the LS formalism (blue short-dashed) considerably departs from the other curves. 

Including higher order terms in equation (|44|l so as to better account for the void-in-cloud effect does not change the 
results for big voids. This is expected from our model in which big voids are unlikely to be embedded in a larger collapsing 
region. However, including the void-in-cloud effect modifies the ratio of void abundances for the smallest voids. Namely, a 
positive fni tend to decreases the number of small voids. This may originate from an increase of high mass halos, which is 
such that small voids are more likely to sit inside an collapsing region. Our prediction is at best qualitative because higher 
order terms are needed to describe accurately the small voids regime (Our results shown in figure [5] indicate that the third 
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Figure 6. Ratio of void abundances of f„i ^ to the /„; = case. Left and right panels show f n i = 100 and —100 respectively. Curve 
labels are the same as the previous figure, with the exception of the solid (black) showing the PS prediction (the square brackets in 
equation H28I I) and an addition dot-long-dashed (orange) curve showing the same equation but ignoring the d(aS3)/ds term. Models 
without the void-in-cloud effect (black solid, blue short-dashed, and orange dot-long-dashed) include the scale dependence of (JS3; the 
others use |ctS 3 | = 0.022. 



order approximation (magenta) is valid only for v > 0.6). We have not pursued the inclusion of higher order corrections 
here since constant barriers may not be a good approximation for the formation of halos and voids as suggested by the 
analysis of (Gaussian) initial conditions of cosmological simulations. Nonetheless, our formalism allows the incorporation of 
scale dependence barriers, which shall be useful when numerical measurements are available. 



5 DISCUSSION 

We extended Doroshkevich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics 
to the local non-Gaussian /„; model. We showed that, up to second order in this is straightforward because, at fixed 
overdensity, the distribution is the same as when /„; = (i.e., Gaussian initial conditions). Our analytic formulae are in good 
agreement with measurements of the distribution of (Ai, A2, A3) in Monte Carlo realizations of the f n i distribution (Figures[T] 
and [2}. 

Our extension of Doroshkevich's formulae to the local non-Gaussian f„i model provides the first step in the study of triaxial 
structure formation in models with primordial non-Gaussianity. This is interesting because, for Gaussian initial conditions, 
halo formation is more triaxial than spherical. In particular, in the triaxial collapse model, the evolution of a patch depends 
on its initial overdensity 8 as well as the parameters e and p, which describe its initial ellipticity and prolateness. We showed 
that the distribution of e and p at fixed S is unchanged from when /„; = (equation [22}. Therefore, equation (|31[1 . which was 
determined for / n j = models, should continue to be useful even when /„; 7^ 0. 

We applied our formulae for the initial shear field to study the change in halo and void abundances in the local non- 
Gaussian model (Section [3]). When /„; = 0, halo abundances predicted by a model in which halo formation is associated with 
having all three initial shear eigenvalues above some critical value are in better agreement with the simulations than those 
implied by the usual overdensity threshold criterion (Figure [3]). However, the predicted dependence of f(y, /«;)// (f, 0) f n i is 
in better agreement with the simulations (see Fig. [4| in the latter case. 

To understand this better, we extended the elli psoidal collapse formalism to the local no n-Gaussian model. We included 



the moving barrier formulation of ellipsoidal collapse (Shet h et al 1 l200ll ; ISheth fc Tormen|[2 002 ) using two different approaches. 



The first is analogous to that used in the case of spherical collapse (Section [A}; the predicted dependence on f„; did not agree 
with m easurements from the simulations. The second is an extension of the excursion set approach following lLam fc Sheth 
(2009a). For the case of a constant barrier (associat ed with spherical collapse), this approach reproduces the results of 
Lo Verde et al. 1 200Sl ) and Maggiore fc Riottol 1 2009b ). Its extension to moving barriers appears promising since it matches 



the measured halo counts when /„; = as well as the dependence on / n ; (Figure [4}. 
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For /„; models, differences in t he density field evol ved from Gaussian and non-Gaussian initial conditions are more 
dramatic in the underdense regions l|Lam fc She th 2009b). Figure [5] shows a number of predictions for void abundances for 
fni — 0, while Figure [6] shows the effect of using our triaxial formalism to study how these depend on /„;. The trends are 
generally like those in the halo abundance, exce pt that the dependen c e on the sign of /„; is reversed (compare Figures [4] 
and[6|. This is consistent with the recent work of Kamionkowski et al. ( 2009h . Still, one might have expected that an excess 
of massive halos also implies an excess of large voids (since all the mass is concentrated in a smaller volume). This is what 
is indeed found for /„; = when one increases the normalisation amplitude as of the fluctuation field. However, this is not 
true in /„( models. In this respect, void abundances provide complementary information to cluster abundances, so a joint 
comparison of both could be used to put constraints on the level of non-Gaussianity. 

We also demonstrated how the void-in-cloud effect can be included in models where / n ; 7^ 0. We used the constant barriers 
as an example. We found that the inclusion of the void-in-cloud effect modifies the abundances of small voids. As a result, 
models with positive f n i show a strong decrement in very small voids whereas models with negative /„; show the opposite. 
This may due to the enhancement of high mass halos for /„; > which effectively increases the probability of finding small 
voids inside a high mass halo. Higher order terms in equation (|44[) and a more accurate description of scale dependent barriers 
will be needed to make more quantitative predictions. 

Our results have other applications which we have not completed. The distribution of the eigenvalues of the initial shear 
field can be used to study the shapes of halos and voids; combining the signals with the halo/ void abundance and shape 
distribution would further constrain the value of /„;. This is a subtle point because, although halo shapes are expected to 
correlate with the parameters e and p, we have shown that, at fixed S, the distribution of e and p does not depend on /„;. 
Hence, naively, the shape distribution is not informative. In practice, one usually averages over a range of halo masses. In 
other words, these will have a range of S/a values, so the result of this averaging may depend on f n i, for the same reason 
that the distribution of e (equation [23} depends on f n i. An d finally, we are in th e process of extending our nonlinear redshift 
space probability distribution function for the dark matter ( Lam fc Sheth 20081 ') to these /„; models. 
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APPENDIX A: CRITICAL VALUE APPROACH 



The derivation of equation (|26[) , which assumes that F is simply related to P, without writing the intermediate steps associated 
with the excursion set approach, is ad hoc. For example, if we assume the same logic that leads to equation (|26[1 but choose 
the integrand p(S\R) so that it returns equation (|29p in the case /„; = 0, then we find that 
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The first line of this expression has not made its way into the /„; = literature - a testament to how much more popular 
the excursion set approach (see below) has become. It does not provide a particularly good description of the measured 
dependence on /„; , so we do not consider it further. 



APPENDIX B: BIVARIATE EDGEWORTH EXPANSION FOR VOID ABUNDANCES 

We used the bivariate Edgeworth expansion to study the first crossing probability when /„; 7^ in Lam fc Sheth ( 2009a! ) . We 
can define the analogy of Q mn here: 
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where the h mn were defined in Lam fc Shethl (|2009al ) 
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with po(x) = e~ x / 2 /\/27r. In addition, we may define 

Q S = [° dSp (8 + 8 v ,s\Sc,S)H 3 (8c/Vs)=H 3 (8c/Vs) 1+eli[{5v zMMF ER. 
If we again ignore the scale dependence of a S3, then 
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where £(s, S) is the same as in lLam fc Shethl l|2009af ) with b = 8 V and B = 5, 



